Two-Cooper-pair problem and the Pauli exclusion principle 
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While the one-Cooper pair problem is now a textbook exercise, the energy of two pairs of electrons 
with opposite spins and zero total momentum has not been derived yet, the exact handling of Pauli 
blocking between bound pairs being not that easy for N — 2 already. The two-Cooper pair problem 
however is quite enlightening to understand the very peculiar role played by the Pauli exclusion 
principle in superconductivity. Pauli blocking is known to drive the change from 1 to A'^ pairs, but 
no precise description of this continuous change has been given so far. Using Richardson's procedure, 
we here prove that Pauli blocking increases the free part of the two-pair ground state energy, but 
decreases the binding part when compared to two isolated pairs - the excitation gap to break a 
pair however increasing from one to two pairs. When extrapolated to the dense BCS regime, the 
decrease of the pair binding while the gap increases strongly indicates that, at odd with common 
belief, the average pair binding energy cannot be of the order of the gap. 

PACS numbers: 74.20.Fg, 03.75.Hh, 67.85. Jk 



I. INTRODUCTION 

The first step towards understanding the microscopic 
grounds of superconductivity was made by Frohlichi who 
has realized that electrons in metals can form bound pairs 
due to their weak interaction with the ion lattice, which 
results in an effective electron-electron attraction. A few 
years later. Cooper has considered^ a simplified quantum 
mechanical problem of two electrons with opposite spins 
and zero total momentum added to a " frozen" Fermi sea, 
i.e., a sea of noninteracting electrons. Within the Cooper 
model, an attractive interaction between these two elec- 
trons is introduced, this interaction being localized in a 
finite- width layer above the "frozen" Fermi sea. Cooper 
has shown that such an attraction, no matter how weak, 
leads to the appearance of a bound state for the two ad- 
ditional electrons. This result was demonstrated for a 
single pair although it was fully clear that conventional 
superconductivity takes place in a macroscopic system of 
electrons paired by such an attraction. 

One year later, Bardeen, Cooper and Schrieffer— (BCS) 
have proposed an approximate solution of the quantum 
many-body problem for electrons with opposite spins at- 
tracting each other. A very important result of the BCS 
theory is the existence of a gap in the excitation spectrum 
above the ground state. In the BCS model, the potential 
layer, in which an attraction between electrons with op- 
posite spins acts, extends symmetrically on both sides of 
the Fermi level. This implies that indeed a macroscopic 
number of electrons interact with each other. In order 
to avoid the difficult problem associated with the Pauli 
exclusion principle between a given number of same spin 
electrons, the grand canonical ensemble was used. The 
original formulation of BCS theory is also based on a 



variational ansatz for the ground state wave function: 
The wave function is taken with all the electrons feeling 
the attraction, paired, i.e., "condensed" into the same 
quantum-mechanical state. 

It was, however, emphasized by Schrieffer that electron 
pairs are not elementary bosons because they are con- 
structed from two elementary fermions^, so that their 
creation and destruction operators do not obey simple 
bosonic commutation relations. Schrieffer also claimed 
that the large overlap which exists between pairs in the 
dense BCS configuration cuts any link with the two-body 
Cooper model, the isolated pair picture thus having little 
meaning in the dense regime^. In spite of this claim, it is 
rather obvious that the many-electron BCS configuration 
can be reached from the one-Cooper pair limit by sim- 
ply adding more and more electron pairs into the layer 
where the attraction acts, until the layer becomes half- 
filled. A canonical procedure of this kind would allow 
one to see the evolution of correlated electron pairs from 
the dilute to the dense regime and to understand deeper 
the role of the Pauli exclusion principle in fermion pair- 
ing. Notice that such an approach can also be considered 
as a useful and well-defined toy model for the crossover 
between local and extended pairs of attracting fermions, 
which in the present time attracts large attention within 
the field of ultracold gases^"— . The crossover problem 
is still open even for the simplest case of the " reduced" 
BCS potential for fcrmion-fcrmion interaction: a varia- 
tional solution has only been proposed long time ago by 
Eagles^i and also by Leggetl^. It also uses a BCS-like 
ansatz for the ground state wave function. 

A possible way to tackle the problem in the canoni- 
cal ensemble, i.e., for a fixed number of electron pairs, is 
to use the procedure developed by RichardsoniSili. It al- 
lows us to write the form of the exact TV-pair eigenstatc of 
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the Schrodinger equation in the case of the so-cahed " re- 
duced" BCS potential which is the simplest formulation 
of the electron-electron interaction mediated by the ion 
motion. The eigenstate, as well the energy of N pairs, 
read in terms of N parameters Rn, which are 

solutions of N nonlinear algebraic equations. Although 
Richardson's approach greatly simplifies the problem by 
avoiding a resolution of a A^-body Schrodinger equation, 
the solution of these equations for Ri, Rn in a com- 
pact form for general N remains an open problem. One 
of the difficulties is due to the fact that N is not a pa- 
rameter in these equations but only enters through the 
number of equations. This is rather unusual and makes 
the A''-dependence of the system energy quite uneasy to 
extract. Nowadays, Richardson's equations are tackled 
numerically for small-size superconducting granules con- 
taining countable numbers of electron pairs^^. We wish 
to add that the canonical approach has also been used 
in the form of a variational fixed-A^ projected BCS-like 
theories, see e.g. Re£i^. 

The goal of the present paper, is to extend the origi- 
nal Cooper's work for one electron pair to two pairs: we 
analytically solve the two Richardson's equations in the 
large sample limit. Our work can be considered as an 
initial step towards the establishing of the precise link 
which exists between dilute and dense regimes of pairs, 
since it indicates a general trend for the evolution of the 
ground state energy with the increase of pair number, 
i.e., overlap between pairs. Richardson's equations are 
here solved by three methods. They of course give iden- 
tical results but shine different light on these equations. 
The approaches to tackle Richardson's equations, pro- 
posed in this paper, in fact constitute perspectives for 
the extension to a larger number of pairs and hopefully 
to the thermodynamical limit. 

The solution we obtain shows that the average pair 
binding energy is smaller in the two-pair configuration 
than for one pair. This result can be physically un- 
derstood by noting that electrons which are paired are 
fermions; therefore, by increasing the number of pairs, we 
decrease the number of states in the potential layer avail- 
able to form these paired states. The energy decrease we 
here find is actually quite general for composite bosons^^. 

However, extrapolation of this understanding to the 
dense BCS configuration faces difficulty within the com- 
mon understanding of BCS results. Indeed, it is gen- 
erally believedi^ii^ that the pair binding energy in the 
dense BCS limit is of the order of the superconducting 
gap A. At the same time, this gap is found as expo- 
nentially larger than the single pair binding energy ob- 
tained by Cooper. According to the tendency we here 
revealed, the average pair binding energy in the dense 
regime should be smaller than that in the one-pair prob- 
lem. 

This discrepancy motivated us to focus on what is 
called pair binding energy and more generally "Cooper 
pair" in the various understanding of the BCS theory. 
Usually, pairs are said to have a binding energy of the 



order of A. However, such pairs are introduced not ab 
initio, but to provide a physical understanding of the 
BCS result for the ground state energyi^iii. Pairs with 
energy of the order of the gap are called "virtual pairs" 
by Schrieffer— . They represent couples of electrons ex- 
cited above the normal Fermi level for noninteracting 
electrons, as a result of the attraction between up and 
down spin electrons. Since the Fermi level is smeared 
out on a scale of A by the attraction, the number of 
such pairs is much smaller than the total number of elec- 
tron pairs feeling the attraction. The latter were named 
"superfluid pairs" by Schrieffer—. By construction, the 
concept of " virtual pair" breaks a possible continuity be- 
tween the dilute and dense regimes of pairs in a somewhat 
artificial way. By contrast, staying within the framework 
of "superfluid pairs" greatly facilitates the physical un- 
derstanding of the role of Pauli blocking in superconduc- 
tivity as well as in the BEC-BCS crossover problem. Our 
results in fact demonstrate the importance of a clear sep- 
aration between the various concepts of "Cooper pair" 
found in the literature. 

We wish to mention that the results presented in this 
paper do not have straightforward experimental appli- 
cations. The main goal of this paper is to reveal the 
general trend for the evolution of energy spectrum when 
changing the number of pairs and to make a first step to- 
wards a fully controllable resolution of the A^-pair prob- 
lem. However, even a two-pair configuration has a rela- 
tion to real materials having correlated pairs of fermions, 
because this configuration corresponds to a dilute regime 
of pairs, realized in some systems. Conceptually, the 
overlap between pairs can be tuned either by changing 
fermion-fermion interaction or total number of pairs. We 
here show that by increasing the overlap between pairs, 
we block more and more states available for the construc- 
tion of paired states. For the first time, a dilute regime of 
pairs was addressed by Eaglcs^i in the context of super- 
conducting semiconductors having a low carrier concen- 
tration. In particular, it was shown in this paper that the 
excitation spectrum in the dilute regime is controlled by 
the binding energy of an isolated pair (in agreement with 
our results) rather than by a more cooperative gap which 
appears, when pairs start to overlap. Thus, this picture 
is quite similar to the isolated-pair model considered by 
Cooper. 

There also is a variety of unconventional superconduc- 
tors which are characterized by rather short coherence 
length that implies pairs not overlapping so strongly as 
in conventional low- Tc materials. For instance, it was 
argued in RefJ^ that the BEC-BCS crossover might be 
relevant for high- Tc cupratcs. Some experiments seem 
to support this idea, for example Refi^ where experi- 
mental data on the dependence of the superconducting 
transition temperature on Fermi temperature are col- 
lected for various superconducting materials. This anal- 
ysis indicates that conventional low- Tc superconductors 
stay apart from short-coherence length materials, includ- 
ing heavy fermion superconductors. Thus, it was argued 
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that to understand these unconventional materials, it is 
appropriate to focus on the most basic aspect, i.e., on the 
short coherence length, rather than to introduce more ex- 
otic and less generic concepts^S. It was also shown that 
the very recently discovered Fe-based pnictides, which 
constitute a new class of high- Tc superconductors, should 
be understood as low-carrier density metals resembling 
undcrdopcd cuprates^i, so that it is possible that the 
BEC-BCS crossover phenomenon is relevant for these 
materials as well. Quite recently, it was demonstrated 
in Refi^ that size quantization in nanowires made of 
conventional superconductors can result in a dramatic 
reduction of the coherence length bringing superconduct- 
ing state to the BEC-BCS crossover regime. We finally 
would like to mention that the two-correlated pair prob- 
lem has received great attention within the ultracold gas 
field, sec e.g. Rcf.— . All these examples demonstrate that, 
paradoxically, the Cooper problem seems to be more rel- 
evant to modern physics than several decades ago. It 
is also worth mentioning that BCS Hamiltonian, which 
only includes interaction between the up and down spin 
electrons with zero pair momentum, is oversimplified. 
Nevertheless, fermionic pairs in the BEC-BCS transition 
regime have not been described yet in a fully controlled 
manner even within this Hamiltonian. One of the possi- 
ble strategies to tackle this crossover therefore is to find 
a precise solution of the problem for the simplest Hamil- 
tonian and only after that, to turn to more elaborate 
Hamiltonians. 

The paper is organized as follows. In Section II, we 
briefly recall the one-Cooper pair problem to settle no- 
tations. In Section HI, we present two solutions to the 
two-pair ground state, as well as a discussion of the pos- 
sible excited states. We conclude in Section IV. In the 
Appendix, we give another exact solution to the two-pairs 
Richardson's equations which shines a different light to 
the problem. 

II. THE ONE-COOPER PAIR PROBLEM 

Let us briefly recall the one-Cooper pair problem. We 
consider a Fermi sea \Fq) made of electrons with up and 
down spins. An attractive potential between electrons 
with opposite spins and opposite momenta acts above 
the Fermi level epo- This potential is taken as constant 
and separable to allow analytical calculations. In terms 
of free pair creation operators = '^ilt'^-ki reads as 

V= -F^?i;k'W7k/3^/3k (1) 

k%k 

y is a positive constant and = 1 for epo < £fc < 
epo + ^■ 

We add a pair of electrons with opposite spins to the 
"frozen" sea |fo)- When the pair has a nonzero momen- 
tum, it is trivial to see that ap-|-a^p/^ |^o) with p ^ p' 
is eigenstate H = + V where Hq = J2k,s ^kols^ks, 



its energy being {sp + Sp'). If the pair has a zero total 
momentum, the H eigenstates are linear combinations of 
Pk\Fo)- We look for them as 

|Vi)=5]G(k)/3t|Fo) (2) 

k 

The Schrodinger equation — f i) = imposes 
G(k) to be such that 

[2ek - £i] G(k) - Vwu ^ «;k'G(k') ^ (3) 

For 2ek ^ £i, the eigenfunction G(k) depends on k as 
Wk/(2£k — £i), so that \tpi) is only made of pairs within 
the potential layer, as physically expected. The eigen- 
values such that 2ek S\ for all k within the potential 
layer then follows from Eq.(3) as 

k - ti 2 ]^^^ 2e - fi 

po is the mean density of states in the potential layer. 
This leads for a weak potential, i.e., a dimensionless pa- 
rameter V = pqV small compared to 1, to 

£i ~ 2£fo ~ £c (5) 

Ec 2^6-^/" (6) 

As seen below, it will be physically enlightening to 
rewrite this one-Cooper pair binding energy as 

e,^{pon)i2e-^^ypo)=Nnev (7) 

Nn = Po^ is the number of empty pair states in the po- 
tential layer f2 from which the Cooper pair bound state 
is constructed, these states being all empty in the one- 
Cooper pair problem. £y = 2e~^/'"/po appears as a bind- 
ing energy unit induced by each of the empty pair states 
in the potential layer, ey only depends on the potential 
amplitude V and the density of states po in the potential 
layer. 

Eq.(7) already shows that the wider the potential layer 
ri, the larger the number of empty states feeling the po- 
tential from which the Cooper pair is made and, ulti- 
mately, the larger the binding energy Sc- We can also 
note that the pair binding energy depends linearly on the 
number of states available to form a bound state. This 
remark is actually crucial to grasp the key role played by 
Pauli blocking in superconductivity: Indeed, this block- 
ing makes the number of empty states available to form a 
bound state decrease when the pair number increases - or 
when one pair is broken as in the case of excited states. 

III. THE TWO-COOPER PAIR PROBLEM 

We now add two pairs having opposite spin electrons 
and zero total momentum to the Fermi sea \Fq) and we 
look for the eigenstates {H — £2) \ip2) = as 

|V2) = EG(ki.k2)/3i,/^ll^o) (8) 



4 



The bosonic character of fermion pairs which leads to 
= /^L-^ki' allows us to enforce G'(ki,k2) = 
G(k2, ki) without any lost of generality. The Schrodinger 
equation fulfilled by G'(ki, k2) is somewhat more compli- 
cated than for one pair. To get it, it is convenient to note 
that 



wpl^l \Fo) 



(9) 



the factor {l — Sk^k^) being necessary for both sides of the 
above equation to cancel for ki = k2. When used into 
{H - £2) |-02) = projected upon (Fo| ^ki^ka we get 



= (1 - 5k,kJ [(2£ki + 2£k.-f2) G(ki,k2) 

-V Wki ^ WkG(k, k2) + (ki o k2) 

V k^^ks 



(10) 



The above equation makes G(ki,ki) undefined. This 
however is unimportant since the ki = k2 contribution 
to IV'2} anyway cancels due to the Pauli exclusion prin- 
ciple. For ki 7^ k2, the equation fulfilled by G(ki,k2) 
follows from the cancellation of the above bracket. With 
probably in mind a ( ki, k2) decoupling , Richardson sug- 
gested to split £2 as 



£2 ^ Ri + R2 



(11) 



with Ri 7^ i?2, a requirement mathematically crucial as 
seen below. We can then note that 



2eki + 2£k2 ^ £2 



1 



1 



(2£ 



i?i) (2£k. - R2) 2eki - ^1 2ek. - R2 



(12) 



with (i?i,i?2) possibly exchanged. This probably led 
Richardson to see that the symmetrical function con- 
structed on the LHS of the above equation, namely 

is an exact solution of the Schrodinger equation provided 
that Ri and R2 arc such that 



1-^E' 



Wk 



2V 



Ri 



Ri — R2 



= (i?i O R2) (14) 



as obtained by inserting Eq.(13) into {H — £2) \'4'2) = 0. 
Note that the denominator in the above equation clearly 
shows why (i?i,i?2) are required to be different. The 
fundamental advantage of Richardson's procedure is to 
replace the resolution of a 2-body Schrodinger equation 
for G(ki,k2) by a problem far simpler, namely, the res- 
olution of two nonlinear algebraic equations. 

This procedure nicely extends to N pairs, the equa- 
tions for Rn reading as Eq.(14), with all possible 
R differences. However, to the best of our knowledge, the 
analytical resolution of these equations for arbitrary TV 
has stayed an open problem, even when N = 2. We now 
show how we can tackle this resolution analytically, first 
through a perturbativc approach, and then through two 
exact procedures. 



A. Perturbative approach 

A simple way to tackle the Richardson's equations an- 
alytically is to note that Eq.(4) allows to replace 1 in the 
LHS of Eq.(14) by the same sum with Ri replaced by 
£1. If we now add and substract the two Richardson's 
equations, we get two equations in which the potential V 
has formally disappeared, namely 



E 



Wk 



2£k — i?l 2£k — i?2 



2E 



2£k - £1 



E 



Wk 



2£k — ^1 2£k — i?2 



Ri — R2 



(15) 



(16) 



V is in fact hidden into £1 . This is a wise way to put the 
singular V dependence of Cooper pairs into the problem, 
at minimum cost. 

In view of Eq.(15), we are led to expand the sums 
appearing in Richardson's equations as 



E 



= E' 



2£k — i?i 2£k — £1 + £1 — Ri 

k k 



E '^n{Ri — £iy^ 



(17) 



n=0 



where Jq — 1/V while J„>o is a positive constant given 
by 



J„ = E 



Wk 



PO In 



(2£k - £1 



j^^— (18) 

^ 2 n£" ^ ^ 



-2n/v 



(19) 



for V small. For this expansion to be valid, we must have 
\Ri — £i\ < 2£k — £1 for all k. This condition is going to 
be fulfilled for large samples, as possible to check in the 
end. 

It is convenient to look for Ri through Ci = 
{Ri - £1) /Ec with i = (1, 2). Eqs.(15, 16) then give 



E-(^r 

n—l 



G? 



(20) 



00 J 

(Gi - G2) J2 -i^i - ^2 ) - -27c (21) 



The above formulation evidences that the Richard- 
son's equations contain a small dimcnsionless parameter, 
namely 



7c = 4/7Vc 



(22) 



where = po^c- Indeed, Nc is just the pair number 
from which pairs start to overlap. This makes Nc large. 
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and consequently 7c small compared to 1, in the large 
sample limit. 

For 7c = , the solution of the above equations reduces 
to Ci = C2 = 0, i.e., 82 = 2fi. The fact that the two-pair 
energy £2 differs from the energy of two single pairs 2£i 
is physically due to Pauli blocking, but mathematically 
comes from a small but nonzero value of 7c. 

To solve Eqs.(20, 21) in the small 7c limit, it is conve- 
nient to set Ci = S + D and C2 = S — D. This allows us 
to rewrite Eqs.(20, 21) as 



lc/2 



[/i + f (D2 + 352) + . ..]+S[l2+hiD^ + S^)- 



pairs. It is in fact possible to avoid this 7c expansion as 
we now show. 

Through the perturbative calculation, we have found 
that the difference Ri — R2 is imaginary at first order in 
7c. It is possible to prove that this difference is imaginary 
at any order in 7c: The Richardson's procedure amounts 
to add an imaginary part to the two-pair energy, in order 
to escape into the complex plane and avoid poles in sums 
like the one of Eq.(4), the two "Richardson's energies" 
then reading a,s Ri = R + iR' and R2 = R — iR\ with R 
and R' real. 

(23) 

R is by construction real since Ri -\- R2 ~ 2R is the 
energy of the two Cooper pairs. In order to show that R' 
. . . J also is real, let us go back to Eq.(16). In terms of {R, R'), 
this equation reads 



+ + Ia(^d^ + qd^s^ + S^) 



(24) 



Their solution at lowest order in 7c reads 7c/2/i 
-1)2 ~ 2Sh/l2. When inserted into Rx + R2 = 2£i 
{Ci + C2)ec, this gives the two-pair energy as 



h 

£2 ~ 2£i -I- 7c —2 Ec 
2If 



2fi + — (l + 2e-2/" 
Pa ^ 



(25) 



Using the expression of £1 given in Eqs.(5,6), we can 
rewrite this energy as 



2£Fo + — ) -£c (l--r^ 
2£Fo + — ) -ev{Nn - 1) 



(26) 



Compared to the energy of two single pairs 2£i = 
2{2ep^ — Ec), we see that Pauli blocking has two quite 
different effects, (i) It first increases the normal part of 
this energy as reasonable since the Fermi level for free 
electrons increases. The first term in Eq.(26) is nothing 
but 2eF„+2{eF„ + -^): one pair has a kinetic energy 2eFg 
while the second pair has a slightly larger kinetic energy 
2(eFn -|- — ), the Fermi level increase when one electron is 
added, being 1/po- (ii) Another less obvious effect of the 
Pauli exclusion principle is to decrease the average pair 
binding energy. Indeed due to Pauli blocking, (TVo — 1) 
pair states only are available to form a bound state in 
the two-pair configuration, while all the Nq pair states 
are available in the case of a single Cooper pair. 



B. Exact approach 

The perturbative approach developed above, through 
the (Ri — £1) expansion of the sum appearing in the 
Richardson's equations, helped us to easily get the ef- 
fect of Pauli blocking on the ground state of two Cooper 



E 



Xl + i?'2 i?'2 



(27) 



where = 2ek — i? is real. We then note that this 
equation also reads 



i?'2| 



= R 



,/*2 



\R' 



i?'2| 



(28) 

From it, we readily see that, since the LHS and the 
bracket are both real, R'*^ must be real. R'*^ can then 
either be positive or negative, e.i., R' can be real or imag- 
inary, which produces (i?i,i?2) either both real or com- 
plex conjugate. 

To show that (i?i, R2) cannot be both real, we go back 
to Eq.(16). By noting that ^ wy^ is nothing but the num- 
ber N^-i of pairs in the potential layer, we can rewrite this 
equation as 



1 



1 



2ek-i?i 2£k-i?2 Nn{Ri-R2) 
Ak 



Ri ^ R2 



E 



(2ek-i?i)(2£k-i?2) 



(29) 
(30) 



where 



Ak = (i?i - R2f + ^(2ek - i?i)(2ek - R2) 



(31) 



It is possible to rewrite the second term of Ay using 4a6 
[a + 6)2 - (a - 6)2. This leads to 



Ay 



(1 - ^)(i?i - R2f + ^(4£k -Ri- R2? (32) 



Since the number of pairs TVq in the potential layer is far 
larger than 1, Ay would be positive if (i?i, i?2) were both 
real. For (i?i,i?2) outside the potential layer over which 
the sum over k is taken, the sum in Eq.(30) would be 
made of terms with a given sign, so that this sum cannot 
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cancels. Consequently, solutions outside the potential 
layer must be complex conjugate whatever 7^. 

For (i?i,i?2) complex conjugate, i.e., R' real, the sum 
over k in Eq.27, performed within a constant density of 
states, leads to 

= I arctan — — — — arctan — —] (33) 

2R' V R' R' ' 

where Xp^ ~ 2eFo ^ ^- If n.O'^ take the tangent of 
the above equation, we find 



The major advantage of this exact procedure is to 
clearly show that the above result corresponds to the 
dominant term in both, the large sample limit by drop- 
ping terms in hi front 1/po, and the small poten- 
tial limit by dropping terms in e~^/^ in front of As 
seen from the first expression of R in Eq.(35), the Pauli 
exclusion principle induces a double correction, in 1/Nc 
and in to the one-pair binding energy ec = 2fle~'^^^ . 
However, the corrections in 1/Nc ends by giving a poten- 
tial free correction to the 2-pair energy £2 because, in a 
non-obvious way, it in fact comes from a simple change 
in the free electron Fermi sea filling, as seen from the 
second expression of R in Eq.(35). 



tan ■ 



2nR' 



PoR' R'^ + Xf,{2VL + Xf„) 



(34) 



Turning to Eq.(15), we find that it reads in terms of 
(E, R') as 



E 



XI + i?'2 



1 

V 



(35) 



If we again perform the integration over k with a constant 
density of states, this equation gives 



R' 



{2n + Xp, 



R 



12 



(36) 



R and R' then appear as the solutions of two algebraic 
equations, namely Eqs.(31) and (33). Unfortunately, 
they do not have compact form solutions. 

It is however possible to solve these equations ana- 
lytically in the large sample limit, po then goes to in- 
finity so that Nq, and Nc are both large. In this limit 
tan(2/poi?') ^ 2/pqR' to lowest order in (1/po)- Eq.(31) 
then gives i?'^ ~ Xi.„(217 + XF^)/Nn. 

For po infinite, i.e., iV^ infinite, R' reduces to zero so 
that, due to Eq.(33), z — XFa/{2Vt + Xpo) reduces to 



-2/v 



Eq.(33) can then be rewritten as 

2 _ -i/v 1 + z/Nn 



1 + l/zNn 



(37) 



Since e "^^^Nn — Nc/2 is also large compared to 1, this 
gives the first order correction in 1/po to , 



e-Vv [1 _^ (g-2/t. _ e^/'')/2N<:i]. From Xf„ = 2nz/{l 
z) which for z small reduces to Xpo — 2Vl{z + z"^), we end 
by dropping terms in e~^/", with R at first order in 1/po 
given by 



R 



2eFo - 2f7e-2/'' ( 1 



Nc 



1 . Is 

Pa Nil 



1 



(38) 



Since E2 = 2R, this result is just the one obtained from 
the pcrturbative approach given in Eq.(26). 



C. Excited state 

We now consider the 2-pair excited states with a bro- 
ken pair having a nonzero total momentum, as possibly 
obtained by photon absorption. Such a pair does not feel 
the BCS potential, so that it stays uncorrelated. These 
excited states thus read 

1^1; k,k') = E^(ki)/?t^a|,^alk,^ |Fo) (39) 

To derive the equation fulfilled by F(ki), it is conve- 
nient to note that, for k 7^ k , 

= (5kip(l-(5k,k-(5kik')a|,^alk'il^o) (40) 

the bracket insuring cancellation for ki = k or k', as 
necessary due to the LHS. It is then easy to show, from 
the Schrodingcr equation [H — Ei^ kk') |V'i; k,k') = pro- 
jected upon (Fol a_k'iakt/^p that 







[(2ep + ek + ek' -fi, kk')^'(p) 



-Vw, 



q#k.k 



u;qF(q) 



(41) 



This makes i^(p) undefined for p = k or k'. This is unim- 
portant since the corresponding contribution in [V'l; k,k') 
cancels due to the Pauli exclusion principle. For p 7^ (k. 
k') the equation fulfilled by i^(p) is obtained by enforc- 
ing the bracket of the above equation to cancel. Follow- 
ing the one-Cooper pair procedure, we get the eigenvalue 
equation for one broken pair (k, — k') plus one Cooper 
pair as 



1 
V 



E 

p#k,k' 



2£r 



- £k + ek' 



^1, 



kk' 



(42) 



A first possibility is to have the two free electrons in the 
two lowest states of the potential layer, namely Ek = Efq 
and Ek' = Sfo + 1/po- The p-state energy in the above 
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equation must then be larger than e^p^ with e^'' = £fo 
n/po; so that Eq.(39) merely gives 



1 

V 



Po 



(2) 



2de 



2e + 2eFo + l/Po - kk' 



(43) 



By writing + as e^^^ + rj(") with fl - n/po, 

Eqs.(4, 5) for the single pair energy readily give 



Si 



kk' 



2eF„ + -)c2 (sF„ + -) 2 (n -) e-y- 
PoJ V PoJ V PoJ 

(44) 

Another possibility is to put the two free electrons in 
the second and third lowest states of the potential layer, 
namely Ek = EFc, + 1/po and ek' = EFc, + 2/po- The p- 
state energy in Eq.(39) can then be equal to epa or larger 

than ep\ In this case, Eq.(39) gives 



1 
V 



1 



Po 



2£p„ — E 



2 /,(3) 

■ 'Pa 



2de 
2e-E 



(45) 



in which we have set = fi, kk' ^ 2eFo ^ S/po- By E as 
2e^j — 2n^^^e~'^^'"x, the above equation gives x through 



Ec (1 — 1/-/Vn) of the average pair binding energy when 
going from one to two correlated pairs. Since Nq = p^fl 
is far larger than Nc ~ Po£c, the gap increase however is 
far larger than the binding energy decrease. 

The changes we obtain in the excitation gap and in 
the average pair binding energy when going from one to 
two pairs, are a strong indication that the gap in the 
dense BCS configuration cannot be simply linked to the 
pair binding energy, as commonly said. Indeed, the pair 
binding energy is going to stay smaller than Ec = 2Q,e~'^l'" 
due to Pauli blocking in the potential layer, while the 
experimental gap in the dense regime is known to be of 
the order of Oe~^/" which is far larger than Ec- 

We wish to stress that, in addition to the excited states 
considered in this section, in which the broken pair ends 
by having a non-zero momentum, there also are excited 
states, not included into the present work. In these ex- 
cited states, the two pairs still have a zero momentum 
but correspond to i?'s located somewhere in the quasi- 
continuum spectrum of the one-electron states, i.e., in- 
between two one-electron levels. For such i?'s, it is not 
possible to straightforwardly replace summation by inte- 
gration in the Richardson's equations as we did through- 
out the present paper. 



xN'!. 



Log 



1 + xe-2/'' 



(46) 



IV. CONCLUSION 



where = 2r2'^^^e~^/"/9o is close to TVc, i-e., large 
compared to 1 in the large sample limit. This gives 
a; ~ 1 -f 21'N ^ , so that the energy kk' would then 
be equal to 



kk' 



7/3 

4eFo + - - 2 O - - 
Po V Po 



(47) 



This energy is larger than the one given in Eq.(41) with 
the broken pair in the two lowest energy levels of the 
potential layer. 

Such a conclusion stays valid for broken pair electrons 
in higher states: the minimum energy for a broken pair 
plus a correlated pair is given by kk' in Eq.(41). The 
excitation gap to break one of the two Cooper pairs into 
two free electrons A = £i. kk' ^ ^2, thus appears to be 



A 



3 

Po 



3 



(48) 



We can then remember that the excitation gap for a 
single pair is equal to EFo + (^Fo + ^ ) ^ (2e Fo ^ Ec) , i-e. , 
£c + The broken pair being again in the two lowest 
states of the potential layer, this brings an additional 
^ contribution to the average pair binding energy Ec- 
Eq.(45) thus shows that the gap increases when going 
from one to two pairs. This increase in fact comes from a 
mere kinetic energy increase induced by Pauli blocking. 
It is worth noting that, while Pauli blocking induces an 
increase of the gap, it produces a decrease from Ec to 



We here extend the well-known one-pair problem, 
solved by Cooper, and consider two correlated pairs of 
electrons added to a Fermi sea of noninteracting elec- 
trons. The Schrodinger equation for the two-pair ground 
state has been reduced by Richardson to a set of two 
coupled algebraic equations. We here give three differ- 
ent methods to solve these two equations analytically in 
the large sample limit, providing a unique result. These 
methods are perspective for the extension to an arbi- 
trary number of pairs in order to hopefully cover the 
crossover between dilute and dense regimes of Cooper 
pairs, as well as to apply them to nanoscopic supercon- 
ductors. Although the two-pair problem we here solve, is 
only a first step toward the resolution of this quite funda- 
mental problem, it already allows us to understand more 
deeply the role of Pauli blocking between electrons from 
which pairs are constructed. We show that this blocking 
leads to a decrease of the average pair binding energy in 
the two-pair system compared to the one-pair configura- 
tion. This decrease is due to the fact that by increasing 
the number of pairs, we decrease the number of available 
states to form bound pairs. 

This two-pair problem actually has some direct rela- 
tion to real physical systems, where correlated pairs are 
more local than in conventional BCS superconductors. 
We can mention underdoped cuprates, heavy fermion su- 
perconductors, pnictides, and ultracold atomic gases. It 
was shown long time ago^, in the context of supercon- 
ducting semiconductors with low carrier density, that the 
excitation spectrum of such dilute system of pairs is con- 
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trolled by the binding energy of an isolated pair rather 
than by a cooperative BCS gap. Hence, this picture is 
very similar to the classical Cooper model. Within the 
two-pair configuration that we here solve using Richard- 
son's procedure, we reach the same conclusion for the 
excitation spectrum. We also reveal how the compos- 
ite nature of correlated pairs affects their binding ener- 
gies through the Pauli exclusion principles for elementary 
fermions from which the pairs arc constructed. 

The extrapolation of the tendency we find to the dense 
BCS regime of pairs, indicates that the average pair bind- 
ing energy in this regime must be smaller than that of an 
isolated pair. At the same time, it is generally believed 
that the pair binding energy in the BCS configuration is 
of the order of the superconducting gap, which is much 
larger than the isolated pair binding energy. To under- 
stand this discrepancy, we must note that there are two 
rather different concepts of " pairs" in the many-particle 
BCS configuration. Those with energy of the order of 
the gap are introduced not ab initio, but enforced to 
have a gap energy in order to provide a qualitative un- 
derstanding for the expression for the ground state en- 
ergy, found within the BCS theory. These entities, called 
" virtual pairs" by Schrieffer— , correspond to pairs of elec- 
trons excited above the Fermi sea of noninteracting elec- 
trons. These virtual pairs have to be contrasted with 
what Schrieffer calls "superfluid pairs"—, made of all the 
electrons with opposite momenta feeling the attracting 
BCS potential, the number of these pairs being much 
larger than the number of " virtual pairs" . Staying within 
the framework of " superfiuid pairs" greatly helps to un- 
derstand the dilute and dense regimes of pairs on the 
same footing. 
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VI. APPENDIX 

In this Appendix, we propose another exact approach 
to Richardson's equations which may turn more conve- 
nient for problems dealing with a pair number larger than 
two. 

We start with Eqs.(14) and calculate the sum by again 
assuming a constant density 



2 _ /-^f-o+f^ 2de 4 



(Al) 



with Ri possibly complex. Instead of Ri, we are going to 
look for a, = (217-f 2eFo -i?j)/(2£_Fo with i = (1,2). 



Since Ri = 2ei?„ — 2ri/(ai — 1). the above equation yields 

2 2 (ai-l)(a2-l) , 

V Nq ai - a2 

where Log denotes the principal value of the complex 
logarithmic function, i.e., the one that satisfies — tt < 
Im (Logai) < TT 

By adding the same equation with (1,2) exchanged, 
we readily get 



aia2 



(A3) 



This equation is nothing but Eq.(33), since for i?i = i?-|- 
ii?', we do have ai = {2Vl+XFa-R-iR')l{XFo-R-iR'), 
with i changed into —i for (i?i, ai) changed into (i?2, 
02). From Eq.(A3), we conclude that ai = e'^/'"t while 



02 



A- 



Next, we note that, since the two-pair energy £2 
Ri + R2; which also reads 



£2 = 4eF„ - 2n- 



ai + 02 



:4/" + l-(ai+a2) 



(A4) 



is real, (ai -f 02) must be real. This implies t + l/t = t* + 
l/t* or equivalently it - t*){tt* - 1) = 0. Consequently, 
t is either real or such that \t\ = 1, i.e., t = e*"^. To 
choose between these two possibilities, we consider the 
difference of the two Richardson's equations as written 
in Eq.(A2). This difference which first appears as 



= Log 



4 (ai-l)(a2-l) 



02 Nq ai - 02 
reads in terms of t as 



-V-t.i-f (t4)Log. 



(A5) 



(A6) 



For t real, {t — t^^)Logt is always positive, except for 
t ~ 1 where it cancels. This shows that the RHS of 
Eq.(A6), equal to 2 for t — 1, stays essentially smaller 
than 2 for far larger than 1. Since e^/" + e^^/^ is 
far larger than 2 for v small, we conclude that Eq.(A6) 
cannot be fulfilled for t real. 

The other possibility is t = e"^ with < |(^| < tt, 
so that Log< — iip We then have ai = 03, i.e., i?i = 
i?2. This shows that the two Richardson's energies are 
complex conjugate, as found by the other exact approach. 
When used into Eq.(A6), this t leads to 



ip sin ip + cos ip = Sc 



(A7) 



-2/1, 



. Once Eq. ( A7) for (p is solved. 



with 6c - ^ 

the two-pair energy given in Eq.(A4) follows from 



£2^4.6 



Fo 



An 



e^/" cos ip — 1 



e^h + 1 - 2e2/" cos 



(A8) 



The solutions of Eq.(A7) cannot be expressed in a com- 
pact form in terms of classical functions. We however 
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see that the two dimensionless terms in Eq.(A7), namely 
2/iVn and (5c, are small. Furthermore is smaller 

than Sc- The function in the LHS of Eq.(A7) is increas- 
ing from 2/Nq up to a maximum ~ 1.82, then decreasing 
down to —2/Nci as ip runs in [0, tt], and still decreasing 
on [tt, 27r]. This shows that Eq.(A7) admits exactly one 
solution in the interval (0, 7r/2), another one in the inter- 
val (7r/2, tt) but no solution in [tt, 27r]. Changing ip in —ip 
would provide also two solutions on (— tt, 0) but this just 
corresponds to exchange Oi and 02; so that they cannot 
be considered as distinct solutions. For these solutions. 
If sin If stays close to zero, so that (p is close to or to tt. 
For (f = 0, the RHS of the above equation reduces, for v 
small, to 4£i?Q — iQe^^^"" which is just twice the energy 
of a single Cooper pair as given in Eq. (6). The effect of 
Pauli blocking on this two-single pair energy results from 
a large but finite number of pair states A^^ in the poten- 
tial layer, as physically expected. We can note that, by 
contrast, ip = tt would lead to E2 close to 4ei?(, -1-406"^/"". 
This solution has to be sorted out because it corresponds 
to i?i and i?2 located in the complex plane very close to 
the real axis where the one-electron levels are positioned, 
so that the distance between them and this real axis is 
of the order of 1/po- This prevents substitution of dis- 
crete summation by integration, in Eq.(Al), as discussed 



above. 

For ip close to zero, Eq.(A7) gives the leading term in 
l/N^ as ~ (gi/" - £-1/^)2/ Afj. The ratio in Eq.(A8) 
then reads for ip small as 



1 



=2/d _ I 



1 - 



^2 e2/.(g2/. ^ 

2 {e^h - 1)2 



1 



^A/v _ _ 1 ^ g-2/. 



(e2/« _ ly 



(A9) 



When inserted into Eq.(A8), we end with 

£2 ^ 4SF„ - 4J7e-2/" + + 2e-^/A (AlO) 

which is nothing but Eq.(26). 

The main advantage of this second exact method is to 
have the two-pair energy reading in terms of ip which fol- 
lows from a single equation, namely Eq.(A7). By con- 
trast, to get £2 through Xpg ~ 2eFa — £2, the 
other exact method, we must solve two coupled equa- 
tions, namely Eqs.(31) and (33). 
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